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We provide a new determination of the charm quark mass using the Highly Improved Staggered 
Quark (HISQ) action, finding mf s (3GeV) = 0.983(23) GeV. Our determination makes exten- 
sive use of second order lattice perturbation theory in matching the bare lattice mass to the MS 
^vq scheme. This matching utilises both traditional diagrammatic perturbation theory and weak cou- 

pling simulations. The second of these techniques allows us to extract perturbative coefficients 
from Monte-Carlo simulations and the process of doing this is laid out in some detail here. 

OO 

o 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 

' Email:iana@triumf.ca 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



Matching mo and 



I. F. Allison 



1. Introduction 

The quark masses are important both as fundamental parameters and, more pragmatically, as 
inputs to experimental determinations of CKM matrix elements [1]. The charm quark is particu- 
larly important because of the large flavour physics program, but it has been somewhat neglected 
on the lattice due to the difficulty in simulating it accurately. In these proceedings, we make use of 
the highly-improved staggered quark (HISQ) action [2] to extract a value for the charm quark mass 
from dynamical lattice QCD. Empirically, the HISQ action is known to reduce the &(a s a 2 ) errors 
which remain in the AsqTad action, and which are thought to be taste-changing errors, making pre- 
cision charm physics possible [2]. It does this by repeating the AsqTad link smearing, which further 
suppresses taste changing interactions, and also by correcting the dispersion relation through ad- 
justment of the Naik term coefficient. In this work, we aim to use these features to calculate the 
charm quark mass from r\ t correlators on the lattice. 

In quoting a determination of a quantity such as the charm quark mass, it is customary to 
convert to the MS renormalization scheme and to use a standard scale (e.g. 3GeV). To do this 
directly from a determination of the bare lattice QCD mass requires lattice perturbation theory. The 
trend towards increasingly complicated actions, such as HISQ, has made most calculations of this 
type a major computational undertaking. One tool developed in response to this is the use of weak 
coupling or high-p simulations [3]. At sufficiently large values of the coupling j3 (equivalently, 
sufficiently small lattice spacings) a lattice simulation will have a small physical volume and a 
very large cutoff(« %jd). These are precisely the conditions required to probe the perturbative 
regime of QCD and when perturbation theory is done in this way (by Monte-Carlo) all orders are 
automatically included. When it is used in combination with a technique like constrained curve 
fitting, the high-/3 technique can allow diagrammatic results to be extended to the next order at 
the cost of running some extra simulations, provided the quality of the high-j8 results is sufficient. 
This approach has been successfully demonstrated in [4]. Of course, there are some complications, 
mostly related to the very small volumes of the simulations, but the two most significant problems: 
the existence of zero modes and Z3 tunneling, are known to be effectively resolved by the use of 
(color) twisted boundary conditions[5], see e.g. [6]. In this work, we will use high-/3 simulations 
with twisted boundary conditions on all of the spatial dimensions to do part of the second order 
matching. 

The recently published determination of [7] also used HISQ quarks and a mixture of contin- 
uum and lattice techniques to calculate the charm quark mass, finding m^(3 GeV) = 0.986( 10) GeV 
The calculation we will present uses a completely different method, extracting m c from r\ c correla- 
tors before manually matching to the MS scheme. Together we view these independent calculations 
as giving important cross checks of one another. 

2. Matching to the MS Scheme 

The lattice charm quark mass am c can be matched to the MS scheme mass m-^ using the 
on-shell mass M as an intermediate stage. 

m M (/i) = M [1 + (B n l + B l0 ) a m (n) + (B 22 l 2 + B 2l l + B 20 ) o^ji)] + ff{o^) 

M = am c [l + {A n L+A w ) a L + (A 22 L 2 +A21L + A20) al] +0(al), (2.1) 
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where L = log am c and / = log \xjM. The relation between mj^ and M (the B coefficients) is given 
to third order in [8] so only the A coefficients are unknown. The connections between ai — > ay, 
and Gtjjs — ► ay are given in [9] and [10] respectively. Writing m-^ in terms of am c to second order 
in Ofy and demanding that the unphysical dependence on L = logam c vanishes gives conditions on 
the coefficients An ,A 22 and A21 which result in the form 



m (n) = am c (1 + (Z n lan +Z W ) a v (aq*) + (Z 22 /^ +Z 2 \l al i +Z20) a v (aq*)) + tf(o#) 

= am c + ci(m ? a)ay(a^*) + (c 2 , ? + C2,g)av(a?*) + ^ , (av), (2.2) 



with Z a „ = log a\i and 



2 _ . 4 15 n/ 

Zll -"7r' Zl °- Al0 "3^' Z22 -2^"3^' 
21ogag 5 \ lllogag 7 2Aio 



221 V 3tt 2 18rc 2 ) n/ % 2 12k 2 k ' 

(logK/aq 4logaq 53 1 \ 

Z2 ° = + + 432^ + 18 J + A2 ° 

2 lllogTT/a^ \ £ 3 2 + log2 22\ogaq 257 

vi,o)Aio + 7^2 wzn ( 2 - 3 ) 



3n 2k l! 7 1U 6tt 2 9 3tt 2 32tt 2 

The splitting of the fermionic and gluonic portions of c 2 in the second line of equation (2.2) is 
motivated by there being only 4 fermionic diagrams for the second order mass renormalization. 
These diagrams have been evaluated using diagrammatic perturbation theory (see [11] for an out- 
line of this calculation, final results are in preparation). The remaining diagrams which contribute 
to C2 tg represent a much larger undertaking and are our motivation for the use of the high-/3 tech- 
nique. One advantage of the split is that we only require quenched results for c 2 . g . The remaining 
unknown coefficients A 10 and A 2 o = A 2 o,g + A20,/, can then be expressed in the following way 

ci 2 

Aio = — + -L, (2.4) 
am c K 

( L 2 l+41og7r ci\ogK/aq \ 7 2 

A20 ' s ~ 6n^ L JnamT) " V ~2^ L + 

2ci ^ 79 + 1321og^-|-24^vi ! o\ 2^C2, g + llcilog^/a^-|- 2^^1,0 
% 12k 2 J 2am r 



A20,/ comes from the diagrammatic analysis of [1 1]. Together with equation (2.3) these coefficients 
allows us to evaluate equation (2.2) and extract a physical values of m-^. 



3. Results 

3.1 The Lattice Charm Quark Mass 

The lattice bare charm quark mass was tuned by adjusting it until the tj c mass agreed with 
experiment on four ensembles of the MILC collaboration's configurations [12]. In this tuning the 
scale was set using MILC values of r\ja and the value r\ = 0.321(5)fm [13] . The bare mass was 
adjusted for any mistuning (a very small effect in all cases) and then converted to the bare tree level 
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mass 1 which is the quantity related to the MS mass in equation (2.2). The m c tree values are given 
in Table 1. 



Size uqchtii uncwi s 


cun c 


'Ic 


m, I GeV 


1 + e 


T\ id 


~ „„ 0.0194 0.0484 
16 3 x 48 

0.0097 0.0484 


0.85 
0.85 


2.26964(17) 
2.27031(16) 


1.100(2) (2) 
1.098(0)(2) 


0.66 
0.66 


2.129(11) 
2.133(11) 


. _ 0.02 0.05 
20 3 x 64 

0.01 0.05 


0.648 
0.66 


1.84153(17) 
1.87142(12) 


1.040(5)(1) 
1.041 (6)(2) 


0.79 
0.79 


2.650(8) 
2.610(12) 


24 3 x 64 0.005 0.05 


0.65 


1.84949(11) 


1.039(3)(2) 


0.79 


2.632(13) 


. M 0.0124 0.031 

28 3 x96 

0.0062 0.031 


0.427 
0.43 


1.30731(11) 
1.31693(12) 


0.9718(5)(11) 
0.9715(5)(11) 


0.885 
0.885 


3.711(13) 
3.684(12) 


48 3 x 144 0.0036 0.018 


0.28 


0.91555(8) 


0.9129(25)(9) 


0.949 


5.277(16) 



Table 1: Simulation parameters for extracting the lattice charm quark mass The value of the correction of the 
naik term, e used here was determined non-perturbatively forcing the "speed of light" to be 1. This differs, 
but not significantly, from the series definition of e used in the perturbative portion of our calculations. The 
value of r\, used to set the scale, was taken to be r\ = 0.321(5) fm 



3.2 High-/3 Perturbative Results 





HISQ 


ASQTAD 


L 3 x T 


6 3 x 16, 8 3 x20, 10 3 x, 12 3 x 20 


6 3 x 16, 8 3 x 20, 10 3 x, 12 3 x20 


P 


15,16,20,24,32,46,62,70,92 


15,16,20,24,32,46,62,70,92 


m 


0.30,0.43,0.50,0.66,0.85 


0.30,0.40,0.50,0.60,0.70 



Table 2: Parameters for the high-/3 simulations. 



We performed high-/3 simulations for valence HISQ and AsqTad quarks at the parameters 
given in table 2. To extract c\ and C2. g we started from the on shell mass M(L,/3) determined by 
simulating quark propagators in a Coulomb+ Axial gauge 2 , and then fitting to the form 

aMp i e = E\ + c\a v {aq*) + c 2 , g ay(aq*) H . (3.1) 

Because high-j8 results are essentially perturbations around the free field, we used constrained 
curve fitting with the first term set to the free field energy of the HISQ action, allowing us to 
evaluate finite volume values for c\ and c%^ g . The values of ay{aq*) were evaluated for each 
simulation by measuring the plaquette and using the three loop expansion of log^ x { given in [4] to 
extract (%v{q* plaq ) which was then evolved to the q* relevant to our simulations. 

3.3 Comparison of HISQ c\ with Diagrammatic Perturbation Theory 

Fits to equation (3.1) were performed including terms up to &{ay ) with priors of ± 5 for all of 
the Cj. We used the resulting c\ values as a check on our method by comparing to the corresponding 

1 The tree level mass is related to the bare mass via m tree = wiq ( 1 — mg + + 5 '^q () — 7 2654400 ra 0° + ® ( m 2 ) 

for HISQ. This relation can be determined from the free field HISQ action. 

2 This is not the traditional maximal-tree gauge but is modified to take account of the twisted boundary conditions. 
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finite volume diagrammatic perturbation theory values as shown in figure 1 . Our results were then 
extrapolated to the infinite volume limit, where we included terms up to fourth order in the fits with 
priors of ± 3 for all parameters. These results were also compared to diagrammatic perturbation 
theory results and are again shown in figure 1 . 



d(L) = d(L = ~) + ^- + ^ + - 



mass 


„MC 
c l 


r PT 
c l 


o 0.50 


0.30 


0.4477(17) 


0.4505(7) 


0.40 


0.43 


0.4932(12) 


0.4921(7) 




0.66 


0.5979(15) 


0.5978(7) 


0.30 


0.85 


0.6634(15) 


0.6693(7) 


0.20 



am=0.30 □ 
am=0.43 
am=0.66 ^-B- 
am=0.85 




Figure 1: c\ infinite volume extrapolations. Boxes are high-j3 results, crosses are diagrammatic PT results. 
The errors come from bootstrapping the entire analysis with 1000 bootstrap re-samples. 



3.4 Comparison of AsqTad C2, g with Diagrammatic Perturbation Theory 

For AsqTad valence quarks, c\ and C2, g have already been calculated using diagrammatic per- 
turbation theory in [14]. We used finite volume values of ci as priors to aid our extraction of C2 lg 
and then extrapolated to the infinite volume limit via 3 

C2, g {L) = c 2 , g (L = oo) + i(X C2A +T C2 ,i logL 2 ) + ^(X C2 , 2 + T C2 , 2 logL 2 ) + • • • (3.2) 

where 7 C2 i = ^X Cl j and X C{ j is the same quantity which appears in extrapolation of c\ and which 
we were able to use as a further constraint in our fits. The results of these fits are given in figure 2 
and we interpret them as lending weight to our C2, g calculation for HISQ. 

4. HISQ c 2 , g From High /3 

The new result which we present here is a determination of the gluonic part of the c 2 coefficient 
for HISQ, for which there are no corresponding diagrammatic perturbation theory results. Again 
we constrained the finite volume c\ coefficients with diagrammatic results and X Cli i from our c\ 
fits. The final results are shown in figure 3. The results are encouraging with the possible exception 
of the result for am = 0.30 which may be affected by finite volume corrections. We are currently 
investigating this possibility by running at larger volumes. 

5. Conclusions 

The high-/3 method has allowed us to extract a second order perturbative coefficient which 
would otherwise have been a very expensive calculation in diagrammatic perturbation theory. To- 

3 The justification for this form for the extrapolation comes from [6]. 
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mass 


r MC 
C 2,x 


C 2,g 


0.30 


1.182(94) 


1.00(2) 


0.40 


1.350(99) 


1.22(3) 


0.60 


1.72(11) 


1.65(3) 


0.70 


2.02(13) 


2.12(4) 



1.40 
1.20 
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Figure 2: AsqTad C2, g infinite volume extrapolations. The boxes are our high-/3 results and the correspond- 
ing infinite volume extrapolation. The crosses are diagrammatic perturbation theory results. For this check, 
the analysis was not bootstrapped, errors are fitting/statistical only. 
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°2,g 


0.30 


0.327(34) 
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0.130(56) 
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Figure 3: C2, g Infinite volume extrapolation for HISQ. The curved lines in the plot are fits to equation 3.3. 
The errors come from bootstrapping the analysis with 1000 bootstrap re-samples. 



gether with calculations of c\, C2, q and equation (2.2) this allows us to provide a two loop determi- 
nation of the charm quark mass using HISQ valence quarks. The final result for each lattice spacing 
we used is given in figure 4. The quoted continuum value comes from fitting all lattice spacings 
simultaneously while demanding a single common charm quark mass, allowing for higher order 
perturbative and discretization errors. 

Beyond statistical/fitting errors, the other important sources of error are the orders excluded 
in the perturbative matching and the overall scale determination. We estimated the error from the 
perturbative matching by repeating our analysis but including the third order perturbative coeffi- 
cients with A30 floated as a very wide prior (0±45). The resulting determinations are also shown 
in figure 4 and suggest that estimating missing order terms by twice a typical value of ay ps 0.22 3 
used in the matching is conservative. The error from setting the overall scale (which we do via n) 
was estimated as 0.5% from an overall error of 1.5% on r\ because the r\ error affects only the 
binding energy for the rj c (see [7]). Including these sources of error, our preliminary result is 

nif{\l = 3GeV) = 0.9830(64) (49) (213) GeV (stat./ntting)(scale)(higher orders). (5.1) 

A more systematic determination, including better estimates of the effect of missing higher order 



6 



Matching mo and m^ 



I. F. Allison 



MSbarmJ3 GeVj/GeV 



mass mf i (3GeV) 



v. coarse 
coarse 
fine 
s.fine 



0.28 
0.43 
0.66 
0.85 



0.9729(53) 
0.9777(12) 
0.9745(59) 
0.9774(18) 



continuum 



0.9830(64) 



1.02 r- 

1.01 - 

1.00 - 

0.99 - 

0.9S - 

0.97 - 

0.96 - 

0.95 - 

0.94 - 



Current-Current analysis 
This Analysis 



0.005 0.010 0.015 0.020 0.025 



Figure 4: Blue points come from full second order analysis, red points come from fits which include a 
parametrization of the third order terms. The black circle comes from the analysis of [7]. 

matching and chiral effects is ongoing and will appear subsequently. At present though, our pre- 
liminary result is in very good agreement with the determination of [7] though with a slightly larger 
eiTor. We interpret this result as a striking demonstration of the capabilities of modern lattice simu- 
lations using highly improved actions such as HISQ to give precise and physically relevant results 
needed by the rest of the particle physics community. 
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